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1 Introduction 


The one-loop QCD vertices have been known for quite some time. In the symmetric case, 
pf = P 2 = P 3 , the one-loop result for the three-gluon vertex in an arbitrary covariant gauge 
was calculated by Celmaster and Gonsalves |1| and confirmed by Pascual and Tarrach |2j] . 
The general off-shell case was considered in the Feynman gauge by Ball and Chiu 
Later, various on-shell results have also been given, by Brandt and Frenkel [j|], and by 
Nowak, Praszalowicz and Slominski || . The most general results, valid for arbitrary values 
of the space-time dimension and the covariant-gauge parameter, have been presented in 
our paper The one-loop quark-gluon vertex (or its Abelian part which is related to 
the QED vertex) was studied by several authors 0 . 

Two-loop corrections to the three-gluon and ghost-gluon vertices have been studied 
in the zero-momentum, limit , which refers to the case when one gluon has vanishing mo¬ 
mentum. In this limit, the renormalized expressions for QCD vertices in the Feynman 
gauge have been presented by Braaten and Leveillc ||. In an arbitrary covariant gauge, 
the relevant results have been presented in our previous paper ||; they will be reviewed 
here. It should be noted that the zero-momentum limit of the three-gluon vertex, as well 
as the relevant limits of the ghost-gluon vertex, are infrared finite, i.e. we do not get any 
singularities of infrared (on-shell) origin. The main argument |9] is just power counting. 

Another limit of interest corresponds to an on-shell configuration, when two external 
gluons are on the mass shell, p\ = p\ = 0. This limit possesses essential infrared (on- 
shcll) singularities The relevant on-shell results for the three-gluon and ghost-gluon 


vertices have been recently presented by two of us 

Information about Green functions is also required for the calculation of certain quan¬ 
tities related to the renormalization group equations, such as the (3 function and anoma¬ 
lous dimensions. The two-loop-order contributions to these quantities were calculated 
by Caswell JT2|, Jones ||R§, Vladimirov and Tarasov fb|, and Egorian and Tarasov [15 


whereas the three-loop-order results were obtained by Tarasov, Vladimirov and Zharkov 
m and by Larin and Vermaseren [17|]. Moreover, recently the four-loop-order expressions 


became available, due to Larin, van Ritbergen and Vermaseren |L8] . 


2 Preliminaries 

The lowest-order gluon propagator in a general covariant gauge is 

( 5 aia2 /p 2 ) (g^ - £ b p^p^/p 2 ) , ( 1 ) 

where = 1 — is the bare gauge parameter. 

The three-gluon vertex is defined as 

= -i 9b f m r wwM (pi,p2,P 3 ), (2) 

where / aia2 “ 3 are the totally antisymmetric colour structures corresponding to the adjoint 
representation of the gauge group (for example, SU(V)); gs is the bare coupling con¬ 
stant. When one of the momenta is zero, the three-gluon vertex contains only two tensor 
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structures ||, 

^A t lM2M3 (Pi Pi 0) (^P/U^P/^ 9Hl^P W 9MUzPM l) -^1 ip ) 

-P^ (W 2 - PmPm/p 2 ) T 2 (p 2 )- (3) 

At the lowest, “zero-loop” order, we have Tf 01 = 1, tA = 0. 

The ghost-gluon vertex can be represented as 

r;r 2a3 (pi,P 2 ;p 3 ) = -i 93 f m pA r fllX 3 (p 1 ,p 2 -p 3 ), (4) 

where p\ is the out-ghost momentum, p 2 is the in-ghost momentum, p 3 and p 3 are the 
momentum and the Lorentz index of the gluon (all momenta are ingoing). For r /i/i3 , the 
following decomposition 0 is useful: 

^ im 3 (P\,P 2 ]P:i) = 9 nM a (P 3 iP 2 iPi) + {terms with p ifl p jfl3 } ■ (5) 

At the “zero-loop” level, f^J] 3 = g MJ/3 , a® = 1. 

In particular, when p 3 = 0 or p 2 = 0 we get 

r hm(~PiP'i °) = 9^M a 3<yP 2 ) + {PmPm 3 term } » 

f w 3 (p, 0; -p) = ^ M a 2 (p 2 ) + {p M P M3 term} , 

as(p 2 ) = a(0,p, -p), a 2 (p 2 ) = a(-p, 0,p). (6) 


The gluon polarization operator and the ghost self energy are defined as 

n w 2 (p) = -^ ai “ 2 (p 2 W 2 -PwPw) Ap 2 ), 


n aia2 (p 2 ) = s aia2 P 2 G(p 2 ) 


-i 


(7) 

( 8 ) 


In the lowest-order approximation J^ 0) = G ( ' 0 ' 1 = 1. 

We shall use dimensional regularization [19|], with the space-time dimension n = 4 — 2e, 
e —> 0. In this paper we adopt the modihcation of the renormalization prescription by 
‘t Hooft [|2(]], corresponding to the so-called MS scheme [^]. In what follows, the notations 
£, a, g (without subscript) correspond to the renormalized (in the MS scheme) quantities. 
In eqs. ([[]), (^|), (£|) they are understood as the bare quantities £#, as, gB- 

The renormalization constants Zp relating the dimensionally-regularized one-particlc- 
irreducible Green functions to the renormalized ones, 


rM (|| 

look in this scheme like 


a,g 


= lim 

£—>0 


Z r ( -,a,9‘ 


r({? 2 }, 


&B i 9 B"> ^ 




= i + £0 

3 = 1 




El 


(9) 


( 10 ) 


where a = 1 — In eq. (|9[) p is the renormalization parameter with the dimension of 
mass. It is assumed that on the r.h.s. of eq. the squared bare charge g\ and the bare 
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gauge parameter as must be substituted in terms of renormalized ones, multiplied by 
appropriate Z factors (cf. eqs. ( 


) and (|I7D ). 

We use the following definitions for renormalization factors: 

= Z-L r MlM 2 /i 3 (pi,P 2 ,P 3 ), 

n£2 oiaa (p) = ^n^(p), 

r (ren) ai a 2 a3 (pi)p2)p3) = ^ f ( Pl , p 2 , p 3 ), 

n( ren ) ai “ 2 (p 2 ) = ^ nai02(p2^ 


The WST identity requires that 


Z 3 /Z 1 — Z 3 /Z\. 


(11) 

( 12 ) 

(13) 

(14) 


(15) 


The results for these renormalization factors in the pure Yang-Mills theory were first 


an ar- 


presented by Jones Jh| (Feynman gauge) and by Vladimirov and Tarasov [|14 
bitrary covariant gauge). The complete results in an arbitrary covariant gauge, includ¬ 
ing the fermionic contributions, were presented by Egorian and Tarasov (cf. also in 
refs. [|,!!). ' _ 

Using (p~5D , the bare coupling constant g 2 B can be chosen (in the MS scheme) as 


9 2 b = [/iV/(4vr )\ g^lZ^Zf = [A7(47r)]y 
where 7 is the Euler constant. The gauge parameter a = 1— f is renormalized as 


as = Z 3 a, so that 
Below we shall also use the notation 
h = g 2 /{^f = a s / ( 47 t ), 


6? = 1 - ^ 3(1 -0- 


where a s = g / ( 47 t) 


(16) 


(17) 


(18) 


3 WST identities 


In a covariant gauge, the Ward-Slavnov-Taylor (WST) identity [^] for the three-gluon 


vertex is of the following form [24 


P 3 3r /4iM2M 3 (Pi>P2,P3) = -J{pl)G(pl) (g,/ 3 pi - p hll pi M3 ) T /13fl2 (p 1 ,p 3 ;p 2 ) 

+J(pI)G(p 2 3 ) (g^ 3 Pl-P2^ P 2 M3 ) f MW (p 2 ,i? 3 ;pi)- (19) 

Consider what follows from (|l^) in the limit when one of the momenta vanishes. 
Contracting with a non-zero momentum, we get 


P™r mMM (p, -p, 0) = -J(p 2 )G(p 2 )a 3 (p 2 ) (g nm p 2 - P„P K ) , 


( 20 ) 


where a 3 {p 2 ) is defined in eq. (| 6 ]). Considering contraction with the vanishing momentum, 
we get a differential WST identity. It can be written in a way which involves just the a 
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functions from the ghost-gluon vertex [0]. For the scalar functions 7\(p 2 ) and T 2 (p 2 ), the 
differential WST identity reads 


Up 2 ) = a 3 (p 2 ) G(p 2 ) J(p 2 ), 


T 2 (p 2 )=2T 1 (p 2 )-2G(0) 


a 2 ( P 2 ) -A (lV 2 J{P 2 ))~P 2 J ( P 2 ) da f P 2 ~ + a 2 ( p 2 ) J{p 

Q' ' nn z 


dp 2 


where the function a 2 (p 2 ) is defined as 


( 21 ) 

( 22 ) 


Up 2 ) = Pu 


d 

9Pl a 


a(P3,-Pi -P3,Pi) 


P1 = ~P3=P 


(23) 


It can be calculated directly at the diagrammatic level (see below). 

Therefore, the differential WST identity makes it possible to define the whole three- 
gluon vertex (not only its longitudinal part) in terms of two-point functions and the 
ghost-gluon vertex. Moreover, it can be used as another independent way, in addition to 
the direct calculation, to obtain results for the three-gluon vertex. 


4 Results for the three-gluon vertex 


The results for unrenormalized one-loop contributions to the scalar functions 7\(p 2 ) and 
T 2 {p 2 ) (in arbitrary space-time dimension) can be found in ref. ||, eqs. (4.30), (4.31), 
(4.33) and (4.34). 

The diagrams contributing to the three-gluon vertex at the two-loop level are shown 
in Fig. 1 of ref. |9|. Note that the non-planar diagrams do not contribute, since their 
colour structures vanish |]25|. When one external momentum vanishes, technically the 


problem reduces to the calculation of two-point two-loop Feynman integrals. To calculate 
the occurring integrals with higher powers of the propagators, the integration-by-parts 
procedure H has been used. For the integrals with numerators, some other known 
algorithms [^, ^7] were employed. Straightforward calculation of the sum of all these 
contributions yields the results for the unrenormalized scalar functions which are presented 
in eqs. (4.8)-(4.11) of ref. 

Using eq. (|IID with 


h r 


Z 1 = l + -[C A (l + !f) - l T ] + h2 [UT 

(-¥-*€+if 2 ) + )( 


+ci 


25 
12e 

45 c 3 c2 


,48 ^32 


(M 

£ - -^-£ : 
^ 16 ’ 


j -c f t 


+ o(h 3 ), 


(24) 


we obtain the renormalized scalar amplitudes appearing in the three-gluon vertex (cf. 

eq- (0)), 


T 1 (ren) = 1 + h [C A (- 


35 


I 20 rp 

9 


+ i£ - i£ : 

18 ^ 2^ 4^ 

4021 _ U _ 2317 f , 15I 113y 2 _Ly 3 I J_y 4 

4^3 f57fi s' 8 S S3 I 144S IfiS ' IfiS 


+h 2 C 2 a (- 

+C a T + 8C 3 + f e - f e) + C f T (f - 16Ca)] + 0(h 3 ), 


(25) 
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ri ren) = h 


Ca (-§ - 2£+ ±£ 2 ) + It] + h 2 [C a T |e 2 ) + 8 C f T 

+ d(-^-Cs + U-^- 287c2 1 19c3 


144 


« -K +o(h 


( 26 ) 

Here and henceforth, we put p 2 = — p 2 in the renormalized expressions. In Feynman 
gauge our expressions agree with those by Braaten and Levcillc ||. 

In eqs. fl25[) and (^6|) we use the standard notation C A for the eigenvalue of the 
quadratic Casimir operator in the adjoint representation, f acd f bcd = C A S ab ( C A = N 
for the SU(iV) group), whereas Cp is the eigenvalue of this operator in the fundamental 
representation (C F = ( N 2 — l)/(2iV) for the SU(iV) group). Furthermore, T = NfT R , 
T r = | Tr(/) = where / is the “unity” in the space of Dirac matrices, Nf is the number 
of quarks, (3 = £(3) is the value of Riemann’s zeta function. 


5 Results for the ghost-gluon vertex 

In order to check the WST identity, we need results for the ghost-gluon vertex in two 
limits corresponding to eq. (||). We shall also need the derivative a- 2 (p 2 ), eq. (|23|) . The 
relevant one-loop results (for an arbitrary n) are listed in Appendix A of ref. |9[]. Two-loop 
contributions to the ghost-gluon vertex are shown in Fig. 2 of ref. |J. Straightforward 
calculation gives the unrenormalized results presented in eqs. (5.6)-(5.13) of ref. 0. 

The derivative (|23|) has been calculated in the following way 0. Let us consider 
the momenta p\ and p% as independent variables, whereas P 2 = —pi — p 3 ■ Therefore, 
the momentum p\ flows from the in-ghost leg to the out-ghost leg. An unambiguous 
Pi path inside the diagram can be chosen as the one coinciding with the ghost line. 
This is convenient, since all we need to differentiate are just two types of objects: ghost 
propagators and ghost-gluon vertices occurring along this path. In this way, we avoid 
differentiating gluon propagators and three-gluon vertices. We also avoid getting third 
powers of propagators. Technically, the propagators and vertices along the ghost path 
were “marked” by introducing an extra argument (say, z). Then, the derivative with 
respect to z was considered, and the rules for differentiating the ghost-gluon vertex and 
the ghost propagator (with subsequent contraction with pi ) were supplied. I 11 this way, 
we just formally differentiate along the ghost line, and then perform all calculations for 
pi = —p 3 = p, P 2 = 0. Finally, extracting the coefficient of g llll3 we arrived at the results 
for the function (^3|) presented in eqs. (5.14), (5.15) of ref. 0. 

Using eq. (|l2|) with 


Zi = l--C A (l-0 + h 2 C 2 A (l-0 


1 (5 — 1A + i / 3 + Id 
c-2 V 8 4 V ^ F V 8 ' 16 V 


+ 0(h s 


(27) 


we obtain the renormalized expressions for the scalar functions occurring in the ghost- 
gluon vertex: 


4 ” n) = 1 + 1 h c A (1 - 0 


+h 2 


ci (f - y 3 - fs - i?c 3 + ie + + \czt] + o(h\ 


(ren) 

a 2 


= l + \hC A £(l-Z) 
+h\ 1-0 


r 1 ' 2 (167 /- 

L A ( ^-fT “8.3- 


43 
144 ^ 


£ - - A£ 

*5 16’ 16 S 


16 v 


-|C^T(l- 5 ) + 0 (h 3 


(28) 


(29) 
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6 Results for the two-point functions 


One-loop results in arbitrary space-time dimension are available e.g. in refs. [|. [28|] (see 
also in Appendix A of ref. 0). Calculating the sum of one-particle irreducible two-loop 
diagrams contributing to the gluon polarization operator and the ghost self energy (shown 
in Fig. 3 and Fig. 4 of ref. [9], respectively), we arrive at the results which are presented 
in eqs. (6.10)-(6.15) of ref. 0. 

Using eqs. 


and (|I4D with 


Z 3 — 1 + 


h r 


£ L 


Ca (| + i{) - IT 


+ h 2 C a T 


l£ 


1 (§-!«)- 5 


2 eJ 


- -C f T 

£ 


+cl 


l(_25_Lltj_ id) i I ( 23 , 15 t _ lti 
f 2 V 12 ' 24 ’ ' 4^ j r \ 8 ' 16^ 8^ 


l£ 


+ o(h 3 ), 


(30) 


z s ~ 1 + -C A (| + \t) + h 2 jc'i — (-1 - + - (§ - 35 $) 


L£‘ 


£ 


we obtain the renormalized expressions for two-point functions 


j(ren) = 1+h 

+h 2 


c ,,(- 

C\(- 


31 


+ £ 


in + f t 


3245 

144 


+ Ca - f C + 2«a + gC 2 - if 3 + g{“) 


(31) 


+C a T (ff + 8<a + f ( - f e) + C f T (f - 160,)] + 0(1, 3 ), (32) 

G<””> = 1 + hC A + h 2 [<7 2 (»2 - gO + |0 2 - gfta) - %C A T] + 0(h 3 ). (33) 

In Feynman gauge, these results agree with those presented in ref. [|J. 

7 Renormalization group quantities 

Using the 1/e term of the renormalization factor Zp (cf. eq. (0)), one can obtain the 
corresponding anomalous dimension yr via 


7r («,r) = -9 2 A C r ] («,/) • 


dg 2 


(34) 


We have checked that in the Feynman gauge £ = 0 (ct = 1) the results for the anomalous 


dimensions 71 , 73 and 73 coincide (in the two-loop approximation) with those from ref. |T 6 
7i = h C A (1+ ) - %T + h* C* A [£+%£- ) - %C A T-4C F T + 


+ 


C A 

(1+ 

f€) 

4 rp 

3 1 

+ h 2 

Cl{ 

H + 

24 ~ 

45 1 
16 ’ 

-r) 

- fC A T- 

AC f T 

IhC 

4(1- 

-0 

+ h 2 C\{ 1 - 

f)( 

-! + 

h) 

+..., 



C A 

(1+ 

h) 

Arp 

3 1 _ 

+ h 2 

<7 ( 

:?+ 

8 ^ 

-\e) 

- 5C a T - 

AC f T 


7i = 

73 = h 

73 = hC A (i + 1{) + h 2 [cl (If - gy - \C a T 


+ .... 
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One can easily check that 71 — 73 = 71 — 73 (this follows from the WST identity ( |H>| ) and 
the definition (0)). Moreover, since f3(g 2 ) = g 2 [271 — 73 — 273 ] (cf. in ref. [[L 6 |) we obtain 
the same result for the two-loop /3 function as those given in refs. I P [TOl l . namely: 
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y(0 


= h 


. 11 - 


c A + |r] + h 2 [- f c\ + f c a t + 4 c f t] + e> (V ). (35) 


Higher terms of the (3 function are available in refs. H16HT8 


8 Conclusion 


We have discussed the calculation of two-loop QCD vertices in an arbitrary covariant 
gauge, mainly in the zero-momentum limit 0. For the three-gluon vertex, we needed to 
calculate two scalar functions, 7\(p 2 ) and T 2 (p 2 ), associated with different tensor struc¬ 
tures, cf. eq. ([|). Two independent ways of calculating these scalar functions have been 
realized. One of them is just the straightforward calculation of all diagrams (cf. Fig. 1 of 
ref. P|). Another way is based on exploiting the differential WST identity. In this way, 
we obtain representations of the scalar functions 7i(p 2 ) and T 2 (p 2 ), eqs. (flTl) and (^), in 
terms of the functions occurring in the ghost-gluon vertex, its derivative (^), the gluon 
polarization operator and the ghost propagator. We have calculated all these functions 
and confirmed the result of the straightforward calculation. 

We have constructed renormalized expressions for all Green functions involved. Note 
that in the zero-momentum limit the three-gluon vertex has no infrared (on-shell) singu¬ 
larities, this is a “pure” case for performing the ultraviolet renormalization. In particular, 
the Z factors (e.g., in the MS scheme) can be constructed just by using the condition of 
finiteness of the renormalized Green functions. 

In principle, the techniques for calculating scalar integrals corresponding to the two- 
loop QCD vertices in the general off-shell case are already available [29]. However, some 
work is still needed to construct the complete algorithm. 
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